Axial Vector Couplings of the Nucleon in Chiral Quark Model 
Incorporating U(1)a Anomaly Effects* 



Xiaoyuan Li 1 ^ and Yi Liao 2 * 
1 Institute of Theoretical Physics, Chinese Academy of Sciences, P. 0. Box 2735 Beijing 100080, 

People 's Republic of China 
2 Department of Modern Applied Physics, Tsinghua University, Beijing 100084, 

People 's Republic of China 



Abstract 

Renormalization of the axial vector currents due to Goldstone loops is studied 
in a simple extension of Manohar - Georgi chiral quark model which incorpo- 
rates U(1)a anomaly effects. The polarized strage quark sea in the polarized 
nucleon results from different renormalization of the flavor singlet and octet 
currents and is in reasonable agreement with the experiment. 
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The chiral quark model of Manohar and Georgi [|IJ offers an explanation of why the 
nonrelativistic quark (NRQ) model works so well for hadrons made up of light quarks. The 
electroweak properties of constituent quarks such as axial vector couplings and magnetic mo- 
ments, however, remain as undetermined coefficients in nonlinear chiral Lagrangian. Using 
the above model in the chiral limit, Weinberg argued some time ago that to the leading 
order in the large N c expansion (N c is the number of colors), constituent quarks behave like 
bare Dirac particles, i.e., with isospin triplet axial vector coupling = 1 and anomalous 
magnetic moment k = 0. The next to leading order corrections have also been estimated by 
sum rules or Feynman diagram calculations in the nonlinear a model or the linear a model 
using the a field as an effective regulator 

What about the flavor singlet axial vector coupling g\ of the quark? Particular interest 
in this problem arises when analyzing the EMC f?| and more recent SMC [|J and SLAC || 
measurements of the nucleon spin structure functions. The original EMC result combined 
with data on nucleon and hyperon f3 decays suggested that quarks carry only a small fraction 
of the nucleon spin and the strange quark sea is substantially polarized in the direction 
opposite to the nucleon spin, thus in conflict with the simple picture of the NRQ model. 
Recent experiment results have reduced the polarization effect of the strange quark sea but 
increased its statistical significance. As emphasized by Kaplan and Manohar [|l(]], upon 
understanding this discrepancy, the key is to remember that the constituent quarks of the 
quark model are not the same things as the current quarks of QCD. Thus, when going from 
the current to constituent quarks the axial vector currents are subject to renormalization 
due to spontaneously broken chiral symmetry in QCD. Furthermore, this renormalization 
in the picture of chiral quark model arises from loops of Goldstone bosons. If the singlet 
and octet axial vector currents are renormalized differently, it is possible to have a polarized 
strange quark sea in the polarized nucleon. In this letter we will study the whole three axial 
vector couplings G^ 8 ' of the nucleon (i.e., the 3rd and 8th components of the octet and the 
singlet) in a simple extension of chiral quark model. We will see that a polarized strange 
quark sea so obtained can indeed be in reasonable agreement with experiment. 



Any theoretical approach trying to give a reasonable value for the flavor singlet coupling 
G° A is faced with the difficulty of how to incorporate the U(1)a anomaly effects. In the sector 
of pure pseudoscalar nonet it has generally been agreed how to do this in a chiral and large 
Nc expansion though controversies still persist [[TTJJfS]]. With the inclusion of quarks, i.e., 
in chiral quark model, simply adding the anomaly term to chiral Lagrangian would double- 
count the anomaly effects because the anomaly is still hidden in the quark integral measure 
when we quantize quark fields. Fortunately there is a way to avoid the dilemma. The key 



point is that one may use a {7(1 ^-neutral quark field |13|]. This is possible because the U(1)a 

transformation of a quark field is just a chiral analog of the ordinary phase redefinition. For 

example, if the quark field if) and the flavor singlet pseudoscalar i] transform under U(1)a 

as if) — > if)' = exp(iwy 5 )ip and i] — > i]' = r/ — uf \/6, one may define the {7(1) ^-neutral 

quark field to be Q = exp(i ^^J 5 )^. The quark integral measure in terms of Q is also 

v6/ 

neutral and no longer generates U(1)a anomaly so that one can now unambiguously include 
the anomaly terms in the chiral lagrangian for pseudoscalars, C ps . For our purpose here 
the relevant point is that the i] Q gets a large mass which is nonzero even in the chiral limit. 
Symmetry considerations dictate the form of chiral Lagrangian at the lowest order in chiral 
and large Nc expansion 

C = Q{ip - M)Q + Q{f + g A Als + g° A £ Q U>)Q 
- QdmiiiP- + | f m|Wp + )Q + £ ps , 

i = exp(z7r a A72/), i = exp^o/Ve/o), 

1 (1) 
P ± = -(l±7 5 ), m = diag(0,0,m s ) 

K = -(id,? - = ^t^o- 

Several remarks are in order. The second term arises from spontaneous chiral symmetry 
breaking, giving quark a constituent mass M. The explicit £77(3) breaking is induced by 
the m term which also gives rise to symmetry breaking terms in £ ps . We have ignored the 
current mass of u and d quarks. The remaining qa and g A terms account for the rule that 



in effective field theory we must include all terms that are allowed by symmetry and of the 
same order by chiral power counting. Compared to the Manohar - Georgi chiral Lagrangian 
the above Lagrangian contains an additional term proportional to g A - The singlet and 
octet axial vector couplings of quark, g A , qa are free parameters in chiral quark model. As 
the Weinberg's argument for g& = 1 in large N c limit was later challenged and there were 
indications that gA = 1 might not be a necessary result of large N c QCD |l4 Hi"6[l , we would 
like to take gA as free and calculate its renormalization effects arising from chiral loops which 
are believed to be one of the important contributions at the next order. In the case of g A 
the situation is more obscure. There is no analogous sum rule to constrain it as used for gA, 
so we set it free as well. 

We are now ready to compute G^ 3 ' 8 , which are defined by the nucleon matrix elements 
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of the QCD current = qj^js—q with A a (a = 1 — 8) the Gell-Mann matrices and 



A° = W~diag(l 1 1) 



< N(p + ?)L7Jb|N(p) >= -±=u(p + q^GAiq 2 ) + q,K(q 2 )} l5 u(p), 

< N(p + q)\jl 5 \N(p) >= ^u(p + q)[ lfl G A {q 2 ) + q^H A (q 2 )] l5 r 3 u(p), 

< N(p + q)\jl 5 \N(p) >= ^u(P + q)h,G s A (q 2 ) + q,H 8 A (q 2 )] l5 u(p), 



(2) 



Note that for G° A we have to specify its renormalization scale \x < A x sb (chiral symmetry 
breaking scale). In the spirit of effective field theory, at scale \x < A x sb we may make the 
appropriate substitution j® 5 — > J^ 5 , where J^ 5 is derived from chiral Lagrangian, 

= lQ^(& a $ - £* a 0Q + 9 -fQlMi^ + 1^1)3 + (meson terms), 

4 o _ 4 (3) 



J% = ^%QlvH^Q + (meson terms). 



V6 

So, effectively we have 



< N\j a jN >=< N\J a jN >=< N|X|N > . (4) 



In the above second equality we have included renormalization effects from chiral loops in 
the coefficients of J^ 5 while the matrix elements themselves are to be evaluated at tree 



level in specific quark models. Since SU (3) is explicitly broken J° 5 and J^ 5 mix under 
renormalization so that a' involves components besides a. Since isospin is conserved, J^ 5 is 
renormalized multiplicatively. We use the NRQ model to evaluate the matrix elements, so 
J^ 5 involves only U and D quarks, 

J 3 J 5 = ^ 2 (Ul,l 5 U-D W D), (5) 

where ao,3,s are effective axial vector couplings of quark in chiral quark model and related 
to observables G A S ' 8 by 
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G A = a , G A = -a 3 , G A = a 8 . (6) 

Explicit calculation of Feynman diagrams in Fig. 1 shows that a 0i 3 5 8 have the following 
structure, 

a = g A (l-A), 

a 3 = g A (l-A)-(2B* + B K ), (7) 
a 8 = g A (l-A)-3B K . 

Instead of writing down the lengthy formulae for A and B^ K \ we emphasize the following 
features. A is a sum of terms contributed by the whole pseudoscalar nonet, while B^ K ^ 
only receives contributions from 7r ± ( K ± , K° and K°). The r/' contributes to ao,3,8 m the 
same way as the Goldstone octet does. This is because we have actually treated rj' as if it 
were a Goldstone boson. Although it is guided by large N c arguments, numerical analysis 
will tell us whether it is a good approximation. The singlet and nonsinglet couplings are 
renormalized differently even in the chiral limit and with g A = g° A . This difference arises 
because Fig. 1(c) and (d) appear only in the nonsinglet channel. Physically it is responsible 
for the polarized strange quark sea in the polarized nucleon, 

AS=«^ = l(g° A -g A )(l-A)+B". (8) 

5 



The splitting between 03 and ag is due to explicit SU(3) breaking, so in the SU(3) limit 
we should have a 3 = a 8 . Indeed, using the 'experiment' value a 3 = 0.75 and a 8 = 0.6 we 
estimate that SU (3) breaking effects are within 30%. 

To quantify our discussion we regularize as usual ultraviolet divergences by the cutoff 
A x sb — 47r/, where / = 84 MeV is the decay constant of the Goldstone bosons in the 
chiral limit. The relevant input is, M = 350 MeV, m s = 150 MeV, m n = 135 MeV, 
m K = 492 MeV, m v = 547 MeV, = 958 MeV, 9 = -20° (77 - 7/ mixing angle), 
f-w = fri = fn 1 — 130/\/2 MeV, fx = 160/v2 MeV. Our discussion does not depend on the 
details of the input. Then G° A 3,8 are functions of qa and g A . As they are sensitive to qa we 
choose a number for gA so that we may get a not-too-bad number for G\ in a reasonable 
range of g A . The result with gA = 1.13 is plotted in Fig. 2. We see that the global pattern 
of G° A 3,8 with a positive g A is in reasonable agreement with experiment. For example, at 
9A = g \ = 1.13, we have G\ = 1.23 ,G 8 A = 0.45 , G° A = 0.16 and AS = -0.10. Considering 
the simplicity of our working Lagrangian this is encouraging. Especially, an important part 
of the strange quark polarization can indeed be attributed to different renormalization of 
the singlet and octet axial vector currents by chiral loops. Although it is possible to fit all 
of G^ 3 ' 8 to experiment, this requires a large negative value for g A . (For example, G A = 0.3, 
G A = 1.25 and G 8 A = 0.6 using g A = 1.16, g\ = —2.14.) A negative sign for the ratio 
( = g A /gA w as also favored by a recent work based on a quantum-mechanical analysis in 



chiral quark model [17]. But we still think this is not very natural because it is hard to 



believe that the next order correction in a good perturbative expansion would change the 
value of G A from —2 or —1 to 0.3 . 

We should mention a weak point in our discussion which seems to deserve further study. 
The corrections to G^ 3 ' 8 computed here are dominated by chiral loops of the Goldstone octet. 
The contribution of the rj' is less important even in the singlet channel while intuitively one 
expects that the 7/ should couple strongly to the singlet channel. This occurs because we 
have actually treat 7/ as a Goldstone boson. We guess that if we put in somehow the 
strong coupling between r( and the singlet current we would have a better expansion of 



the experimental values of G^ 3 ' 8 . This is not totally impossible. An additional diagram 
like Fig. 1(c) or (d) would produce some " B " term which partially cancels the " A " 
term in Oq. Indeed a similar cancellation does occur in the octet coupling 03,8. If the guess 
is really correct we can start with a smaller g A but end up with larger G A ( hence a more 
reliable perturbative treatment in the singlet channel) and G A ' 8 , thus in closer agreement with 
experiment. Inversely, this may imply that some higher order terms in large Nq expansion of 
chiral Lagrangian are probably important. The other point not touched upon in this letter 
concerns the running property of the singlet axial vector current in the region fi < A x sb- 
According to the analysis in Ref. [IJ strong interactions in this region are much weakened, 
compared to the naive extrapolation from the perturbative region of QCD. If we consider 
our computed G A to be evaluated at fi « M and mimic the running in the intermediate 
region A x sb > \i > M by simply using a " scaled-down " QCD ( just as we mimic technicolor 
by using a " scaled-up " QCD ) we estimate G° A (fi = A xSB ) « 0.98 G° A {fi = M) JTj|. 

Finally we argue that g A appearing in the chiral Lagrangian depends in some sense on 
the number of light flavors involved. Usually we work with two flavors ( u and d ) when 
we determine from g A the coupling G\ as measured in the neutron f3 decay. To determine 
couplings of other components, G° A and G A , we surely have to include the strange quark. 
But we should obtain the same G\ whether we work with two or three flavors. To the 
leading order, the relation between g A and G\ is unchanged, e.g. as in the NRQ model. 
Its next order corrections arising from chiral loops are basically determined by quadratic 
Casimir operators of the flavor symmetry group in the symmetry limit. Then one way out 
is that g A also depends on the number of flavors involved. This can be understood in the 
language of effective field theory. When we are working with two flavors we have already 
integrated out the strange quark, K and 77 fields and inserted their renormalization effects 
directly into g A which is a parameter in chiral Lagrangian with two flavors. This explains 
why we required a larger value of g A than usual to fit G\. 

We thank Y. P. Kuang for helpful discussions. 
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Figure Caption 

Fig. 1 Feynman diagrams contributing to renormalization of axial vector currents. Thick 
lines and thin lines represent quarks and pseudoscalars respectively. Solid circles represent 
insertion of current. 

Fig. 2 The computed couplings 3/5G\, G\ and G° A ( upper, middle and lower curves) 
of the nucleon are shown as functions of g A at g\ = 1.13. 
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(c) (d) 
Fig. 1 Xiaoyuan Li, Phys. Lett. B 
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